Renormalization group study of Luttinger liquids with boundaries 
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We use Wilsons weak coupling "momentum" shell renormalization group method to show that 
two-particle interaction terms commonly neglected in bosonization of one-dimensional correlated 
electron systems with open boundaries are indeed irrelevant in the renormalization group sense. 
Our study provides a more solid ground for many investigations of Luttinger liquids with open 
boundaries. 
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Introduction — The electron-electron interaction has 
strong effects on the low-energy physics of (quasi) one- 
dimensional (Id) metals. Such systems cannot be de- 
scribed within Landaus Fermi liquid (FL) theory but 
rather form a different "universality" class: the Luttinger 
liquids (LLs)fi Already decades ago the importance of 
local single-particle inhomogeneities in LLs allowing for 
electron backscattering with momentum transfer 2kp 
(were kp is the Fermi momentum) was emphasized 
Their relevance became even more apparent after it was 
shown that on low energy scales a LL (with repulsive two- 
particle interaction) with a single local impurity in many 
respects behaves as if the chain was cut in two at the 
position of the impurity with open boundary conditions 
at the endpoints^ 5 - This led to a number of studies on 
the physics of LLs with open boundaries&i&^iiii mainly 
using the method of bosonization.ii£ 

The Tomonaga-Luttinger model^ is a translationally 
invariant (no inhomogeneities) continuum model within 
the LL "universality" class for which all correlation func- 
tions can be computed exactly by bosonization. In the 
context of bulk LL physics it plays a role similar to 
the noninteracting electron gas for FL physics.— In this 
model the 2kp two-particle scattering process of spin- 
up and spin-down electrons is neglected. All the other 
low-energy scattering process can be written as quadratic 
forms in the bosonic densities of the right- (around kp) 
and left-moving (around —kp) fermions. Furthermore, 
after linearization of the single-particle dispersion, also 
the kinetic energy is quadratic in the bosonic densitiesi^ 
and the remaining Hamiltonian is that of noninteract- 
ing bosons. Neglecting the above 2kp scattering pro- 
cess leads to an additional conservation law, namely 
the "local" (in momentum space) spin around the two 
Fermi points ikp, which is central to the exact solu- 
tion of the Tomonaga-Luttinger model. In the litera- 
ture one finds two arguments why this so-called gi t ±- 
term can be neglected. Either one follows the original 
idea of Tomonagai^ and assumes a two-particle interac- 
tion which is sufficiently smooth and long-ragend in real- 
space (suppressed screening) such that its Fourier com- 
ponent at 2k p vanishes or, more generally, one adopts the 
so called g-ology renormalization group (RG) approach 
of SolyomJ^ In this it is shown that in an important 



part of the parameter space gi.± is irrelevant in the RG 
sense. The solution gi,± = is a stable fixed point of 
the RG flow.— To (qualitatively) understand the low- 
energy physics the gi^-term can thus be neglected. It 
only affects the numerical values of the other fixed-point 
couplings. 15 

For open boundary conditions the single-particle quan- 
tum number k in which the noninteracting problem is 
diagonal (see below) no longer corresponds to the mo- 
mentum. Therefore the two-particle scattering terms ap- 
pearing in this natural basis have a form different from 
those of the translationally invariant case. Also for LLs 
with open boundary conditions the Hamiltonian contains 
two-particle scattering terms which cannot be written as 
quadratic forms in bosonic densities, when considering a 
general, spin-conserving two-particle interaction. It was 
shown 10 that these terms vanish if Tomonagas rather 
specific assumption of an interaction which is smooth 
in real-space is made. Then open-boundary bosoniza- 
tion can be used to compute all correlation functions 
of the open-boundary analog of the Tomonaga-Luttinger 
model— iIi^iiSi i 10 ! 11 Surprisingly, an RG analysis for more 
general two-particle interactions similar to Solyoms ap- 
proach was so far not discussed in the standard liter- 
ature on open-boundary bosonization. In the present 
Brief Report we close this gap. Using Wilsons "momen- 
tum" shell RG in weak couplin g 16 i 17 we show that all 
two-particle scattering terms naturally arising in a low- 
energy analysis of a Id system with general two-particle 
interaction and open boundaries which cannot be written 
as quadratic forms in the bosonic densities are RG irrele- 
vant (open boundary g-ology analysis). We demonstrate 
that the RG flow equations are the same as the ones of 
the translationally invariant g-ology model. In analogy 
to this the RG irrelevant coupling constants only affect 
the fixed-point couplings of the terms quadratic in the 
bosons and can thus be neglected for a qualitative un- 
derstanding of the low-energy physics of LLs with open 
boundaries. Our analysis puts many of the studies of LLs 
with open boundaries on a more solid ground. 

Model — To be specific we consider the Id electron gas 
on a line between x = and x — L with a general, 
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sufficiently regular two-particle interaction 

V = \ Y, [dx [ dx'i>i(x)i>l,(x') 

xV s<s ,(x - x')ip s ,(x')ip s (x) 

where s, s' denote the spin. The noninteracting one- 
particle eigenstates are given by standing waves (open 
boundaries) (p n (x) — y/2/Lsm (k n x), with quantum 
numbers k n = mr/L, n G N. We emphasize that for 
open boundaries only right-moving electrons with a pos- 
itive Fermi velocity v F = k F /m e appear, where m e de- 
notes the electron mass. The field operator is given by 
ipi^ (x) = y^L-t (p n (x)an}s with the creation and anni- 
hilation operators a^l of electrons in the noninteracting 
eigenstates. We now closely follow Ref. [l^ to rewrite the 
interacting part of the Hamiltonian as 

V = \H H V n° m ,n>, m > a li,s a ln,s> a rn',s>a n >, s 
s,s' n,rn : n' ,m' 

with the matrix elements 

V n,m,n' ,m' ~ i^s^' (k n ~ k n > k rn ~ k m ') 

F s _ s ' (fen k"n' •, k m &m') ~~ -^s.s' (fen ^~ &n') &m ^m') 

-\-F SjS f (k n -(- k n > j k m -\- k ra > 

)]/L, (1) 

where 

F s>s t(q,q') = — [ dx f dx' cos (qx)V StS '(x — x') cos (q'x') . 
L Jo Jo 

(2) 

The kinetic part of the Hamiltonian H = T + V is given 
by 

s n 

with the single-particle dispersion = fc 2 /(2m e ). 

Perturbative Wilson RG — To analyze the Hamiltonian 
H using a Wilson RG we consider the action S(^, 
appearing in the imaginary time functional integral rep- 
resentation of the grand canonical partition function 

= So + S! 

= -^^^2^ s (k n ,co j )[g (k n ,Lj j )}- 1 ^ s (k n ,uj j ) 

s n j 

+f 3 ~ 1 J2 5 i+i,i'+3' v n*m,n',m> 
s,s' n,m,n' ,m' ,j f 

X * s (k„ , LOi ) #y (km , Wj)^ s > (k m > , OJj'^s (k„> ,^i>) , 

with Grafimann variables Matsubara frequencies 

u)j, the noninteracting propagator Qo(k,uS) — [iui — 
f-k + m] (/-* is the chemical potential), and the in- 
verse temperature (3. Here we will only be interested 
in the behavior at zero temperature (/? — > oo) with 



H = ep = k F /(2m e ). In this case we obtain — /i = 
v F (k - k F ) + (k- k F ) 2 /(2m e ). 

We now follow the standard steps of the perturba- 
tive Wilson fc-shell RG;^£ (i) separating the Grafimann 
fields * = + into fast modes \J/> (relative to the 
Fermi point) with A/6 < \k — k F \ < A and slow ones 
(ii) integrating out the fast modes , I /> perturba- 
tively, (iii) rescaling the quantum numbers k = k — k F 
(k' = bk) and the Matsubara frequencies (ui' = 6a;) 18 
and (iv) rescaling the Grafimann fields ^'(k + k F ,uj') — 
^- 1 ^ < (k / /b + k F ,uj'/b). 

We first consider the noninteracting case. Then the 
action can naturally be written as So = S Q < (^> < , \E' < ) + 
S Q > (^' > , 4 r> ) and the second step of the RG procedure 
only leads to a constant in the partition function. There- 
fore the rescaling is the only relevant step. Choosing 
£ = 6 1//2 the iuj — v F (k — k F ) part of the noninteract- 
ing action remains invariant, while the quadratic term of 
the dispersion vanishes as 1/6 for large 6. It is RG irrele- 
vant for the low-energy physics and from now on we work 
with Qo(k,u>) — [iui — v F (k — kp)] . The standard lin- 
earization can thus be justified by the RG. Note that to 
avoid proliferation of symbols we do not introduce a new 
one for the free propagator of the model with linearized 
dispersion. 

Next we study the flow of the quartic part of the action. 
As usual in a weak coupling RG to compute (c /o>) 

where the index refers to taking the expectation value 
with respect to the fast mode part of the noninteracting 
action, we use a cumulant expansion 

He- S ')s> = (Sx)s> + \ (<S?) S > - (Si)l>) +■■■■ 

(3) 

On tree level (first term of the above expansion) and for 
low energy scales (6 — > oo) the RG steps lead similarly 
to the translational invariant case to a two-particle inter- 
action which is purely local in real space V SlS r(x — x') 
V StS '(0)S(x — x'). All correction terms are subleading.^ 7 
Then the integrals in Eq. ([2]) can be performed leading 
to 

F S!S ,(q,q') = V s y(Q) (8 q , q , + 6 q ,- q ,) /2 . 

This has to be contrasted to the case with periodic 
boundary conditions in which only the second Kronecker 
delta appears and the interaction matrix element only 
has the contribution coming from F s>s / (k n —k n i, k m — k m i) 
with k being the momentum^ Since for open boundaries 
all the incoming and outgoing quantum numbers must be 
close to the single Fermi point only the Kronecker deltas 
$n — n',m— m' 3 &n— n',— m+m' 5 and <5 n -i-n'.m+m' in the expres- 
sion Eq. ([T]) for the interaction matrix element contribute. 
As in the translationally invariant g-ology model we now 
generalize the interaction and assume different coupling 
constants for each of the remaining low-energy scattering 
processes and each of the relative spin orientations of the 
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FIG. 1: Topology of the second order diagrams contributing 
to the RG flow. 



two particles. Doing so one can expect that the following 
low-energy RG analysis is applicable not only to the elec- 
tron gas model but also of relevance for a larger class of 
models, including lattice ones. For the interaction matrix 
element Eq. (fl]) this leads to 
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+(\\92S s ,s' + ±92Ss,-s')S n -n',m-m' 
+ ( ||54 S SiS > + ±9i 5 Si - s >^ 5 n -ri ,-m+rn' 



The initial values of the generalized coupling constants 
are functions of the microscopic parameters of the under- 
lying model. The relation between these two parameter 
sets can (in principle) be determined by an additional 
RG step or simpler by perturbation theory. For the cou- 
pling constants we use the same notation-^ gi , and g^ 
as in the translationally invariant g-ology model^ while 
we still make clear that we are dealing with a different 
situation by putting the spin orientation index to the left 
of the ff's3 Two reasons for this become clear already at 
this stage of the discussion. In analogy to the standard 
g-ology the xffi-term is the one which cannot be written 
as a quadratic form in the bosonic density operators of 
the right-moving electrons (for a detailed discussion on 
this, see Ref. llOl ). Also in strict analogy the \\92~ and 
l|5i-terms describe the same scattering processes (as can 
be seen by reordering Grafimann fields and renaming in- 
dices in the expression for S) and could be combined in a 
single coupling constant 1 1 <7 2 — \\9\ (see also the RG flow 
equations below). 

With these interaction matrix elements one now has 
to compute all (second order) connected diagrams of the 
two-particle vertex related to the second term of Eq. ([3|) 
with all the external quantum numbers k set to fc^ and 
the Matsubara frequencies set to zero. Writing down 
the corresponding analytical expressions it becomes im- 
mediately clear that only the one-particle irreducible di- 
agrams give nonvanishing contributions. Furthermore, 
in complete analogy to the g-ology RG of the transla- 
tionally invariant model to leading order the (74-terms do 
not contribute to the RG flow of the coupling constants 
(they do not lead to ln(6)-terms) and do not flow them- 



sclf. They only lead to a renormalization of the Fermi 
velocity, which is a higher order effect if the RG flow of 
the couplings is considered. We thus neglect the 54-term. 
In a straight forward but tedious calculation we evaluated 
all the remaining second order diagrams with the topolo- 
gies as shown in Fig. [TJ The wiggled line stands for any 
of the four coupling constants n7i, n7 2 , and j_72- 

We identified if a certain diagram gives a finite RG flow 
of the vertex and if so to which of the four scattering 
channels it contributes. To perform the remaining fc-sum 
we took the limit L — > 00 (resulting in a semi-infinite 
chain). Setting b = e s and resorting to infinitesimal RG 
steps this leads to the RG flow equations 
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-GrV+iV) , 

-2_l7i (_l7 2 -||7 2 + ||7i) 



(4) 



with 
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for v = ||, _L and I — 1,2. The coupling constants 7 
are now understood to be functions of the infrared cutoff 
parameter s. 

Remarkably, these equations are exactly the lowest- 
order RG flow equations for the coupling constants (di- 
vided by 2irvp instead of nvp) of the translationally in- 
variant g-ology modeljii This provides the complete jus- 
tification for the decision to use the same notation. For 
the solution of the RG equations and an analysis of the 
fixed points we can proceed as it is well documented in 
the literature.— We first introduce the coupling constant 

71,2 = ||7i - ||7 2 + _l7 2 , 
which leads to the RG flow equations 



d.s 
d 

ds 



71,2 



_7i 



-2_l7i 7i,- 



From this it becomes apparent that 7i 2 — J_7i 2 is invari- 
ant under the RG flow which implies that the RG trajec- 
tories form hyperbolas in the ±7i-7i,2-plane. The equa- 
tions can easily be solved and trajectories are sketched 
in Fig. O If initially (for the initial cutoff) all coupling 
constants are small and 'yj 2 > | ^7i J | holds they all stay 
small (the use of perturbative RG is justified) and the 
flow is towards a line of stable fixed points with 



±7i* = 0, 7? 
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2 f-U7if 
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i7 2 * = ±l 2 l + 2 (7i* 2 " 7*,a) 



(5) 
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FIG. 2: Sketch of the trajectories of the weak coupling RG 
flow in the j_7i-7i,2-plane (X — x7i, Y — 71,2)- The dashed 
line indicates a line of stable fixed points. The fixed points 
on the dashed-dotted line are unstable. 

indicated by the dashed line in Fig. [21 Under the above 
restriction on the initial coupling constants the term 
which cannot be written as a quadratic form in the 
bosonic density is RG irrelevant. The initial value j_7i* 
of the j_7i-term only affects the fixed-point values of the 
other couplings pi* The use of the open boundary ana- 
log of the Tomonaga-Luttinger model is then justified 
for studies on the low-energy physics and standard open- 
boundary bosonization&l&^iiii can be used. This allows 
to compute all correlation functionsi^&SJil In the re- 
sulting expressions the coupling constants j_72, 11 72 and 
l|7i must be replaced by the fixed-point values Eq. ([5|). 
For a model in which the coupling constants on the ini- 



tial scale of the weak coupling RG do not depend on the 
relative orientation of the spins of the two scattering elec- 
trons (like the electron gas model) one finds 71,2 = ||7i 
and _l7i = m7i< Therefore the above condition on the 
initial couplings for reaching a stable fixed point reduces 
to the simple requirement that the interaction must be 
repulsive. In this case the trajectory in Fig. [2] flows to 
the origin. 

Summary — We have shown that under quite general 
assumptions on the two-particle interaction the terms 
which are usually ignored in bosonization studies on LLs 
with open boundaries are indeed RG irrelevant. These 
terms which cannot be written as quadratic forms in 
the bosonic densities only affect the fixed-point values 
of the other couplings but do not modify the low-energy 
physics. For this analysis we have used a perturba- 
tive Wilson RG scheme. Although the scattering terms 
for open boundaries are of different nature than in the 
translationally invariant standard g-ology model we have 
shown that the RG equations have exactly the same form. 
Our result puts a large number of bosonization studies 
on LL with open boundaries on a more solid ground. The 
present discussion is limited to the weak coupling regime. 
Numerical studies show that a similar low-energy physics 
can also be found at larger (repulsive) couplings.— We fi- 
nally note that exactly the same results can be obtained 
using the functional renormalization group^i instead of 
the Wilson 
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